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ABSTRACT 

We use atomic hydrogen (H l) data from the Southern Galactic Plane Survey 
to study the kinematics of the fourth quadrant of the Milky Way. By measuring 
the terminal velocity as a function of longitude throughout the fourth Galactic 
quadrant we have derived the most densely sampled rotation curve available for 
the Milky Way between 3 < i? < 8 kpc. We determine a new joint rotation curve 
fit for the first and fourth quadrants, which can be used for kinematic distances 
interior to the Solar circle. From our data we place new limits on the peak to 
peak variation of streaming motions in the fourth quadrant to be ~ 10 km s~^. 
We show that the shape of the average H I profile beyond the terminal velocity is 
consistent with gas of three velocity dispersions, a cold component with Av = 6.3 
km s~^, a warmer component with Av = 12.3 km s~^ and a fast component with 
Av = 25.9 km s~^. Examining the widths with Galactic radius we find that 
the narrowest two components show little variation with radius and their small 
scale fluctuations track each other very well, suggesting that they share the same 
cloud-to- cloud motions. The width of the widest component is constant until 
R < 4 kpc, where it increases sharply. 
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1. Introduction 

From the first detection of the atomic hydrogen (H l) spectral fine in 1951, the A21-cm 
fine has been used as a probe of the kinematics of the interstellar medium (ISM) in the 
Milky Way and the rotation of the Galaxy. The rotation curve of the Milky Way has been a 
long-standing topic of study with important implications for models of the structure, mass 
distribution and dynamics of the Galaxy. In one of the earliest studies of the rotation curve, 
Kwee, MuUer, & Westerhout (1954) noticed irregular depressions in the rotation curve of the 
first Galactic quadrant, which they attributed to a lack of H I between spiral arms. We now 
understand that H I emission does not give a good indication of gas density and therefore the 
dips do not refiect density effects, but instead the kinematics of the H I line provide a direct 
probe of the dynamical effects of spiral structure. In particular, the spiral arms produce 
deviations, both positive and negative, away from circular rotation that are responsible for 
the irregularities observed in the rotation curves of the inner Galaxy (Kerr 1962; Shane 
& Bieger-Smith 1966). These departures from circular rotation have been reasonably well 
explained by both linear and non-linear spiral density wave theories (Lin et al. 1969; Roberts 
1969). 

Interior to the Solar circle, every line of sight passes a single subcentral point where it 
is closest to the Galactic Center and where its projection on the Galactic plane is tangent 
to a circular orbit around the Galactic Center. At the subcentral point the measured radial 
velocity of the H I line reaches its most extreme value, called the terminal velocity. In the 
first Galactic quadrant the terminal velocity is the most positive velocity in H I spectra, 
whereas in the fourth Galactic quadrant it is the most negative velocity. By measuring the 
terminal velocity along all lines of sight between Galactic longitudes —90° and +90° it is 
possible to derive a measurement of the rotation curve of the inner Galaxy along the "locus 
of subcentral points" . 

The shape of the H I velocity profile beyond the terminal velocity can tell us about 
the random motions of H I in the Galaxy. Gas at these velocities is forbidden by Galactic 
rotation and therefore gives us information about the distribution of H I cloud velocities 
at the subcentral point. Kulkarni Sz Fich (1985) employed this technique in an analysis to 
confirm the existence of two populations of H I clouds, those with a velocity dispersion on the 
order of 7-9 km s~^ as found by Burton (1976) and Liszt (1983) and a second population with 
a much larger velocity dispersion of ~ 35 km s~^. The Kulkarni & Fich (1985) observations 
towards I — 30° of a population of "fast" clouds confirmed earlier work by Radhakrishnan 
& Srinivasan (1980) and Anantharamaiah et al. (1984) that had originally identified this 
population of intermediate velocity H I clouds towards the Galactic Center and toward H ll 
regions, respectively. 
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This is the first paper in a series about the kinematics of the inner Milky Way. In this 
paper we develop a new technique for measuring the H I terminal velocity. We apply this 
technique to new, high resolution H I data from the Southern Galactic Plane Survey (SGPS; 
McClure-Grifhths et al. 2005) to determine the H I terminal velocity along the locus of 
subcentral points in the fourth quadrant of the Milky Way. The SGPS data have an angular 
resolution of 130", which far exceeds the resolution of any previous surveys used to explore 
H I kinematics at the terminal velocity. Prom the run of terminal velocities we determine 
a rotation curve for the inner Galaxy. We also use these data to explore the kinematics of 
the neutral hydrogen near the subcentral point. The organization of this paper is as follows: 
in §3 we describe the technique we have used with the SGPS to determine the terminal 
velocity curve and discuss the advantages of using a high resolution survey for terminal 
velocity analysis. In §4.3 we use the SGPS terminal velocity curve to derive the rotation 
curve for the fourth quadrant and we compare these data with a published rotation curve 
of the first quadrant. In §4.3.1 wc derive a new fit to the rotation curve of the inner Galaxy 
incorporating data from both the first and fourth Galactic quadrants and in §4.3.2 we discuss 
the magnitude of the velocity residuals and their relationship to spiral structure. Finally, in 
§5 we discuss the shape of the H I velocity profile beyond the terminal velocity and touch 
on its implications for the kinematics of diffuse H I. 



2. Data 

The H I data discussed here were obtained as part of the Southern Galactic Plane 
Survey (SGPS; McClure-Griffiths et al. 2005). The SGPS is an H I spectral fine and 21-cm 
continuum survey of the Galactic plane between longitudes 253° < I < 358° (SGPS I) and 
5° < I < 20° (SGPS II) and with latitude coverage of |6| < 1?4. The observations and 
analysis of these data are described in full detail in McClure-Griffiths et al. (2005). The 
SGPS data were obtained with the Australia Telescope Compact Array (ATCA) and the 
Parkcs radio telescope and were combined to provide images with sensitivity to all angular 
scales from 130" to several degrees. The SGPS 1 H I data have been released as a scries of 
overlapping sub-cubes, each covering 11° of longitude and 2?89 of latitude with an angular 
resolution of 130" and a spectral resolution of 0.82 km s~^. The noise in the line free channels 
of the H I cubes is 1.6 K. 
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3. Measuring the Terminal Velocity 

On each line of sight near the Galactic plane circular rotation gives a velocity relative 
to the Local Standard of Rest (LSR) that projects on the line of sight as 

vlsr = -Ro sin I {uj - ujq) cos 6, (1) 

where uo is the angular velocity of Galactic rotation at the point of interest, and uoq is the 
corresponding Galactic rotation of the solar circle {R = Rq). Throughout this paper we 
adopt the lAU standard values for the distance to the Galactic Center, Rq = 8.5 kpc, and 
loq — Qq/Ro — 220 km s~^/8.5 kpc. The terminal velocity is defined as the extreme velocity 
on any line of sight at 6 ~ 0°, which occurs at the subcentral or tangent point where the 
galactocentric radius is R — Rosinl for any realistic rotation curve. The projected radial 
component of this maximum allowed circular rotation velocity gives the rotation curve of 
the Galaxy, uj{R), or alternatively Q{R). At this position the helio-centric distance, r, is 
determined from geometry simply as r = Rq cos /. 

For cold gas in a perfectly circular orbit, and in the absence of turbulent motions, the 
H I spectrum would drop to zero at velocities beyond the terminal velocity (in the fourth 
Galactic quadrant: at more negative velocities). However, H I and other real spectral line 
tracers do not drop discontinuously to zero for velocities beyond the terminal velocity because 
the gas has a measurable thermal velocity width as well as bulk random motions. For most 
of the gas in the Milky Way these random motions are roughly 5% or less of the circular 
velocity (Burton 1976), so line profiles near latitude zero drop quite sharply near the terminal 
velocity. Analysis of the shape of the H I velocity profile beyond the terminal velocity can 
give us information about the random motions of the gas in the region near the subcentral 
point. Similarly, understanding the distribution of random velocities helps refine the method 
of estimating the extreme allowed circular rotation velocity from the line profiles. Here we 
carefully fit the terminal velocity profile to reveal both. 

A variety of methods have been employed to define a point along the velocity profile 
as the terminal velocity. For example, Kerr (1962) defined the terminal velocity, Vlsr, as 
the last peak of the profile, whereas Sinha (1978) used a relative definition in which the 
terminal velocity is where the line profile crosses a threshold set at half the value of the most 
extreme allowed velocity emission peak. Shane & Biegcr-Smith (1966) and Celnik, Rohlfs, & 
Braunsfurth (1979) derived functional forms for the shape of the H I profile and the terminal 
velocity and fit their model profiles with a defined terminal velocity to the observed profiles 
to define the exact drop off. More recently, Malhotra (1995) used an absolute brightness 
temperature threshold to define the terminal velocity. We have employed a compromise 
technique, which starts with an absolute brightness temperature threshold to adjust the 
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velocity profile followed with function fitting to define the exact velocity of the underlying 
discontinuity, independent of threshold choice. We made two passes through the data in an 
effort to avoid the effects of an arbitrary choice of measurement parameters. 

3.1. Continuum Blanking 

An important advantage of a high angular resolution survey hke the SGPS is that 
the effects of absorption of continuum emission from beyond the subcentral point can be 
removed from the data. When observed toward a bright continuum source, the H I spectrum 
shows a mixture of emission and absorption at all velocities. The absorption can overwhelm 
the emission, particularly near the subcentral point where velocity crowding causes a long 
segment of the line of sight to pile up in the spectrum near the terminal velocity. For 
a low resolution survey, such as one observed with a single dish telescope, absorbed and 
non-absorbed pixels are averaged with a weighting given by the beam shape. 

The blended emission and absorption in a low resolution survey leads to an underesti- 
mate of the maximum velocity, as absorption towards continuum sources in the beam will 
tend to flatten the steep "cliff edge" of the emission profile near the terminal velocity. A clear 
example of this effect is shown in Figure 1, which shows two spectra at the same longitude 
at two different latitudes (b=±0?2, 1 — 305?25) where the first has bright continuum in the 
background and the second does not. The effect of the absorption is to reduce the height 
of the emission, in extreme cases the brightness temperature can even go negative because 
the continuum has been subtracted from the SGPS data. This causes the apparent rotation 
velocity, whether defined by a threshold or by profile fitting, to be biased to unrealistically 
low values wherever there is bright background continuum emission. To avoid this bias in 
our analysis we have used the SGPS continuum images (Haverkorn et al. 2006) to blank all 
pixels with Tcont > 0.5 K. Overall we blank about 10% of the survey area. Without blanking 
we find that the measured terminal velocities show more scatter from one longitude to the 
next. Single dish surveys are not able to do this because their much larger beam size irrevo- 
cably blends the spectra toward continuum sources with the surrounding emission spectra, 
reducing the measured terminal velocity. 

3.2. Thresholding on Latitude Averages 

After continuum blanking our first step is to chose a brightness temperature threshold 
that selects a velocity that is part way down the steepest slope in the emission profile, the 
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"cliff" just past the subcentral point velocity. An individual spectrum is rather noisy (ctt 
= 1.6 K) so it is necessary to perform some kind of latitude average before thresholding. A 
simple average over latitude would return the mean of this distribution. Alternatively we 
could work with the median, i.e. the 50th percentile T5, or some other value of Tf, defined 
from this distribution such as the 25th or 75th percentile. Because all latitudes are included 
in constructing the spectrum any of these statistics could be used on the distribution of Tj, 
at each velocity to construct a spectrum that has relatively low noise. 

The SGPS interferometer data latitude coverage of ±1° corresponds to 55 times the 
beam width, or about 70 independent points for each beam width in longitude. The scale 
height, h, of the H I gas in the inner Galaxy is about ±150 pc in z (Dickey & Lockman 
1990), which gives ±1° angular width at the nominal galactocentric distance Rq = 8.5 kpc. 
Thus all lines of sight with |6| < 0?5 pass the subcentral point when they are still well inside 
the H I layer (|2;| <^ h), so all spectra at the same longitude show roughly the same sharp 
drop beyond the subcentral point velocity. Corrections to the projected rotation velocity to 
account for the non-zero latitude, both due to the reduced projection of the circular rotation 
velocity vector along the line of sight and due to departures from cylindrical rotation of 
the gas above or below the mid-plane at the subcentral point, are insignificant (~ 10^^) for 
latitudes |6| < 1°. For the terminal velocity analysis described here we have used only spectra 
within |6| < 0?5, so the effect of the finite scale height is very shght. We can therefore consider 
any set of spectra at the same longitude as realizations from the same statistical ensemble, 
differing in their profile shapes due to small scale random motions near the subcentral point, 
but all having the same underlying rotation velocity. 

For each longitude we have used the many independent spectra at different latitudes to 
study the distribution function of H I emission brightness temperature, 3],, at each velocity. 
An example of one such distribution function of at — 0?5 < 6 < -|-0?5, I — 305?25 and 
V = —61.83 km is shown in Figure 2. This is the velocity where the mean of the 
distribution crosses the 10 K threshold, i.e. the most negative velocity with mean Tb{v) > 10 
K. The positions of all four statistics, mean, median, 25th and 75th percentile are marked on 
Figure 2. We have used all four statistics to determine the velocity of threshold crossing. We 
have carried the velocities found with each statistic through the second part of the analysis, 
as discussed in §3.3, to estabfish which statistic is most robust. 

Whatever the choice of threshold in Tf, and the statistic chosen to determine the average 
Tii{l,v) from Ti,{l,b,v), after continuum blanking and averaging we have a function Vt{l), 
which is the velocity of threshold crossing. This is not the rotation curve, but it is the first 
step toward the rotation curve. Using this function, we shifted all spectra in the survey 
by the velocity corresponding to their longitude, Vt{l), and then computed an l-v diagram 
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based on the average in latitude. By shifting the spectra in velocity by Vt{l) we removed the 
gross effects of Galactic rotation so that we could average spectra in longitude. Figure 3 is 
a terminal velocity profile averaged over the longitude range 272° < / < 340°. The velocity 
scale of Figure 3 reflects the shift so that km s~^ corresponds to the threshold crossing, 
Vt{l). The subcentral point is at a positive velocity somewhat offset from vt- The value of the 
offset depends on the threshold choice, and the choice of statistic (mean, median, or other 
percentile) used for the latitude averaging. In the next section we describe how we determine 
the offset to the subcentral point by fitting complementary error functions, erfc{v), to the 
terminal velocity profile. 



The second step in determining the terminal velocity was to fit the shifted spectral 
averages with a functional form to establish a robust velocity for the subcentral point. For 
cold gas with no velocity dispersion the H I profile would drop sharply at the terminal velocity. 
However, the H I profile does not show the ideal edge of a Heaviside step function, but is 
instead smoothed out by the velocity dispersion of the gas. A commonly used technique to 
find the edge of a step function blended by some natural dispersion is to fit a complementary 
error function, erfc, which returns the position of the edge and as well as the width of the 
blending process. The relationship between the step function and the error function can be 
seen clearly if we consider the Heaviside step function as the limit of a complementary error 
function 



where Vo is the velocity offset of the edge and Av characterizes the deviation from an ideal 
edge. The use of an error function to fit the H I profile beyond the terminal velocity is not 
without precedent, Kulkarni & Fich (1985) used this form in their exploration of the velocity 
dispersion of H I. In this application, the width A^;, relates to the Gaussian dispersion as 
Av — -\/2(7„. 

There is also some physical justification for using an error function to fit the H I line 
profile. Beyond the terminal velocity, vlsr, we expect that the gas must become optically 
thin, meaning that the brightness temperature becomes proportional to the column density 
of atoms in each velocity interval, 5v. If the density of gas at allowed velocities is a constant, 
Uo, then the brightness temperature at a velocity beyond the terminal velocity {v < Vo, where 
Vq < for the fourth quadrant) for which the gas is optically thin is 



3.3. Error Function Fitting of the Terminal Velocity Shape 




(2) 




(3) 
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where the upper limit on the velocity integral is not really infinity, but the range of allowed 
velocities is so wide on most lines of sight that the upper limit does not matter. The 
function, f{v — v'), is the velocity distribution function of H I atoms contributing to the 
emission, assumed to be a Gaussian. The constant Co is the familiar 1.823 x 10^^ cm^^ 
(K km s~^)~^ conversion factor between brightness temperature and column density for the 
A21-cm line. Finally, \dv'/dr\ is the magnitude of the inverse of the velocity gradient along 
the line of sight, evaluated at position r that has radial velocity v'. As we can see. Equation 3 
bears a strong resemblance to the complementary error function but for the \dv/dr\ term 
inside the integral. Celnik et al. (1979) showed that for a smooth rotation curve the velocity 
gradient, dv/dr, near the terminal velocity gradient is small and directly proportional to 
distance, r with a slope given by the local value of Oort's constant over R. At the subcentral 
point the velocity gradient crosses zero, such that the absolute value of the velocity gradient 
has a sharp discontinuity at the subcentral point. For optical depths r ^ 1, this sharp 
discontinuity should result in a sharp peak in the H I profile at the terminal velocity. In fact, 
we generally do not observe a sharp peak, but a smoothly rolled off velocity profile, which 
suggests saturation in the line profile near the subcentral point. The effects of saturation 
and spatial variations in the density appear to mask any correlation between the magnitude 
of the velocity gradient and the amplitude of the emission at the terminal velocity, such that 
the erfc function scales to connect smoothly with Tb{vo), whatever it is. Thus if we make 
the assumption that for low longitudes the term Ug \dv'/dv\~'^ is roughly a constant near the 
subcentral point then it can be taken out of the integral, leaving the complementary error 
function for the velocity dependence of the line profile beyond the terminal velocity. Fitting 
this function reveals both the terminal velocity and the width of the underlying Gaussian 
distribution of velocities. 

As was pointed out by Kulkarni & Fich (1985), a single Gaussian for f{v) is not a good 
fit to the data. A single Gaussian is unable to to fit the sharp drop plus a long tail of emission 
beyond the terminal velocity. The next most simple functional description of Ti,{v) is the 
sum of two or three complementary error functions, corresponding to two or three Gaussians 
of widths, (7^: 

3 

T,{v) = J2 a«erfc(-l * {v - Vo)/Av,), (4) 

i=l 

where the width of the error functions, Avi, is related to the original Gaussian width as 
Avi — V^ay^. The top and bottom panels of Figure 3 show fits of two and three error 
functions, respectively, to the longitude average of the shifted Tb{v). Clearly the low-level 
emission tail cannot be well fitted by even the sum of two error functions, instead it is 
necessary to use a sum of three error functions. We discuss the implications of this for H I 
kinematics in §5 below. 
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To determine the true terminal velocity at the subcentral point vlsr = ft + fo, we fit 
error functions to shifted and longitude averaged Ti,{y) profiles to return values for Vo{l)- 
The formal fits were done by minimizing the chi-squared for the fit using Powell's method of 
minimization (Press et al. 1992). We found that fo(/) determined from fitting a sum of two 
error functions was identical to within 2% to Vo{l) determined from fitting a sum of three 
error functions. Therefore, for speed and robustness against noisy spectra we determined 
Vo{l) by fitting the sum of only two error functions. 

The fitted values for Vo{l) depend on the threshold value and the statistic used to 
average the spectra in latitude, but the true terminal velocity, vlsr, is largely independent 
of threshold and statistic. This is illustrated in Figures 4 and 5. Figure 4 shows the threshold 
crossing velocities, Vt, before and after correction by Vo for twenty longitude bins between 

1 — 340° and I — 292°. The four crosses in each longitude bin are the values of Vt obtained 
from the 10 K threshold crossing for each of these four statistics. In each case the outhers are 
the 25th and 75th percentiles, whereas the median and mean statistics give similar results. 
Fitting the velocity offset, Vo, to the corresponding spectrum averages gives different values 
that usually make up for the difference in the percentiles, so that the sum, vlsr = Vt + Vq 
comes out the same for all four statistics. The dotted curve corresponds to v^sR for the 
75th percentile, the dashed curve for the 25th percentile, and the two solid curves for the 
mean and median. For longitudes I < 325° the scatter among the four values of Vt is 1 to 

2 km s""^, but the scatter among the four values of vlsr is only 0.26 km s"^. In this range, 
the continuous curves on Figure 4 are indistinguishable. For longitudes 327?5 < I < 330° 
the points spread apart, but the values of vlsr have less scatter than vt- At all longitudes 
the values of Vlsr. obtained from the mean and median in latitude agree to better than 
0.6 km s~^, except at / = 328?25 where their difference is 4.6 km s~^ and at / = 338?25 
where the difference is 1.2 km s~^. The striking agreement of the values of v^sR. obtained 
from the different choices of statistics to define Vt after correction for the offset velocity, Vq, 
demonstrates that the erfc fitting effectively removes any bias introduced by the arbitrary 
choice of statistic for the latitude averaging. 

The velocity offset correction also effectively corrects for any bias introduced by the 
choice of threshold values. This is shown in Figure 5 where we compare Vt and vlsri both 
derived using longitude averages but with brightness thresholds of 10 K and 20 K. The 
20 K threshold reduces the velocity shifts by 3-5 km s~^, but this is made up by smaller 
offsets between the threshold crossing velocity and the edge of the underlying step in the 
brightness distribution derived from the error function fitting, so that the discrepancy in 
Vlsr is much smaller. The difference in Vlsr. derived using the two thresholds is less than 1 
km s~^ for longitudes / < 325°, and less than 2.5 km s~^ for longitudes / > 332?5. Overall the 
adjustment from the erfc fitting does a good job of correcting for bias due to the arbitrary 
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choice of threshold. The difference in threshold crossing velocity, Vt, is typically 5 to 15 
km s"^ before applying the shift, Vg, obtained from the erfc fits. 

As seen in Eq. 3 and discussed above, the brightness temperature is related to product 
of the local H l density, and the inverse of the velocity gradient, \dv/dr\~^. Because of 
this relationship there is a danger with a pure threshold method that a localised change in 
Uo or dv/dr could increase the local brightness temperature and hence alter the measured 
terminal velocity. The velocity offset correction method that we employ provides some 
robustness against masking density or velocity gradient changes as terminal velocity changes 
by compensating for variations in the threshold velocity. Our estimates of the offset velocity 
made by fitting error functions to the terminal velocity profile are dominated by the velocity, 
rather than the brightness temperature, of the peak of the error function at ~ vlsr- 
Therefore, the value of the threshold, or even its position along the velocity tail, is largely 
irrelevant to our determination of vlsr and we are immune to density and velocity gradient 
effects on Tg. It is worth noting however, that the offset fitting method does not provide 
much protection against very large-scale (factor of ~ 2 or more) variations in Uq or dv/dr 
away from the subcentral point which can effectively shift the velocity of the peak back 
from the terminal velocity. At present there is no method that can resolve these from true 
variations in the rotational velocity and it remains a limitation of all Milky Way rotation 
curves. 

Thus the final values of Vlsr should be repeatable and independent of the choice of 
threshold or statistic to within 0.5 km s~^ for / < 325°, and within 3 km s~^ for / > 332?5. 
In the range 326?5 < / < 332?5 there is a discrepancy in the results from the two thresholds 
that is as large as 10 km s""-*^. Similarly, the range 327?5 < / < 330° gives a large variation 
in threshold crossing velocity using the same threshold but different statistics in latitude 
(see Figure 4). For example, the mean and median for the 10 K threshold give Vt — —143.1 
km s~^ and Vt — —138.3 km s~^, respectively. The offsets derived from the error function 
fitting in that case are Vg — 24.5 km s~^ and Vo = 23.7 km s~^, which reduce the discrepancy, 
but only slightly. This appears to be due to abnormal variation in the slope of the profiles 
beyond the terminal velocity in this region, making the method of latitude averaging, and 
the choice of threshold, give different results. This is possibly evidence for heightened rate of 
stirring up random motions in the medium, due to increased star formation activity in the 
galactocentric radius range 3.9 < i? < 4.5 kpc. As discussed below, this is a region where 
the rotation velocity increases rapidly with R, which may be due to non-circular motions 
associated with a spiral arm. 
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4. The Terminal Velocity Curve 

Using the threshold method and error function fitting we have determined the terminal 
velocity, vlsr-, ^oi 1000 longitude bins of width 3' between / = 339?5 and / = 272°. We used 
a threshold value of 10 K to determine Vt{l) and we fit the sum of two error functions for 
each longitude bin to determine, Vgil). Table 1 gives the measured terminal velocity, vlsri 
and central longitude, Z, for each longitude bin. The terminal velocity curve, together with 
the values for Voil), is shown in Figure 6. For clarity we have not plotted the error bars, 
which should reflect the estimation of the error from on thresholding and fltting described 
in §3.3. 

The longitude range used for this terminal velocity curve analysis has been chosen to 
avoid extreme non-circular motions associated with the inner 3 kpc of the Galaxy. In order 

to translate the measured values of vlsr, to a rotation curve we must assume circular rotation 
along the locus of subcentral points. We know, however, that circular rotation breaks down 
in most of the very inner Galaxy. At < 3 kpc, gas orbits become elliptical because of 
the dynamical effects of the three kiloparsec arm and the bar (Burton & Liszt 1993; Ferriere 
et al. 2007). The dynamics of the H I in the innermost region of the Galaxy will be the 
subject of a future paper and will not be discussed further here. 

Automated function fitting produced spurious fits for ~ 2% of the longitude bins. Some 
of these are visible as cusps in the terminal velocity curve at / = 306°, / = 312° and 
I = 320° in Fig. 6. Examining the H I datacubes we found that at these three longitudes 
there are discrete H I clouds 10-20 km s~^ beyond the terminal velocity, which disrupt the 
fltting. These clouds are spatially only marginally resolved, have low brightness temperature 
features {Ti, ~ 10 K) and FWHM velocity widths of ~ 10 km s^^. They may be similar to 
the outlying clouds described by Stil et al. (2006). In the remainder of the SGPS H I proflles 
used in this analysis the high velocity linewings are smooth and do not show small, discrete 
components. To remove the velocity spikes at / = 306°, 312° and 320° from the subsequent 
analysis we have run a median fllter with a width of 9 samples across the Vo{l) values. 

4.1. Comparison with CO terminal velocities 

H I because of its ubiquity is often hailed as an excellent tracer of the terminal velocity 
curve. However, its ubiquity together with large (~ 5 — 10 km s^^) thermal and turbulent 
line widths can lead to line blending. Because of these limitations, the terminal velocity 
curve has often been studied from observations of the ^^CO J = 1 — transition. Unlike 
H I, the thermal linewidth of the ^^CO line is typically very small, < 1 km s~^ and even the 
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turbulent width is typically only ~ 2 km (Burton 1976), which reduces the problems of 
line blending. Because of this ^^CO determined terminal velocity curves should demonstrate 
less scatter than their H I counterparts. Comparison of ^^CO and H I terminal velocities can 
test how free our H I terminal velocities arc from systematic effects caused by turbulence 
and line-blending. In Figure 7 we have overlaid values for the recent ^^CO terminal velocity 
curve measured by Luna et al. (2006) on our H I terminal velocities as a function of sin/. 
The Luna et al. (2006) values are sampled every 0?125 and were selected as the half-intensity 
point of the extreme velocity peak. 

The agreement between the H I and -"^^CO datascts is striking; the two curves traces each 
other extremely well including many of the dips and wiggles. On average the scatter in the 
^^CO data points is about three times that of the H I data points. The larger scatter in the 
CO data may be partially attributed to the fact that these data are less sensitive than the H I 
data and also that Luna et al. (2006) used a threshold method for determining the terminal 
velocity. However, it is likely that physical effects relating to the small volume filling factor 
of CO contribute significantly to the difference in scatter between H I and CO. The low filling 
factor of CO means that observations at the subcentral point typically detect a very small 
number of discrete CO clouds. The CO observations are therefore sensitive to the individual 
velocity dispersions of the CO clouds, rather than an ensemble average as measured with 
H I. This can lead to a larger scatter for individual measurements. The comparison of these 
two curves raises confidence that H I is not overly dominated by blending or turbulence and 
may in fact indicate that H I is a more robust measure of the terminal velocity. 



4.2. Determining ARq 

Oort's A constant is an important quantity for probing differential Galactic rotation in 
the vicinity of the Sun. Although H I data are not able to determine A directly, terminal 
velocities have been used to measure the product ARq, where 

VLSR^-2ARo{l-\sml\)+v'. (5) 

measured near to the Sun (e.g Gunn et al. 1979). In principle we could determine ARq 
by simply fitting Equation 5 to the SGPS terminal velocities as a function sinZ plotted in 
Figure 7. This is best done near the Sun because A is a, locally determined quantity. However, 
this is complicated enormously by streaming motions both near the Sun and throughout the 
inner Galaxy, which lead us to fit over large ranges of sinZ in order to overcome localized 
non-circular motions. Lin, Yuan, & Roberts (1978) use nonlinear density wave theory to 
demonstrate that the value obtained for ARq depends strongly on the location where it is 



-13- 



measured. Prom Figure 7 it is clear that although the general trend for the SGPS terminal 
velocity is linear, there arc significant departures from a linear fit. We found that the values 
determined for ARq depend on the range of longitudes used in the fitting. We immediately 
eliminated data in the range sinZ > 0.95 because the magnitude of the measured terminal 
velocities at these longitudes (~ 10 km s~^) is comparable to the width of the terminal 
velocity profile as well as any expected deviations due to streaming motions. In Figure 7 we 
show two fits to the data over different ranges of sin I. We found that fits to 0.8 < sin Z < 0.95 
return ARq — 110 ± 2 km s~^, which agrees very well with Gunn et al. (1979) and implies 
A — 13 km s~^ kpc~^ if Rq — 8.5 kpc. However, fits over a broader range of sin I extending 
up to sinZ = 0.35 return ARq = 89 ± 4 km s~^, which imphes A — 10.5 km s~^ kpc~^ if 
Rq = 8.5 kpc. Alvarez, May, & Bronfman (1990) found a similar difference between the slope 
for fits to data up to / = 320° and for data up to / = 340° and attributed the difference to 
anomalous velocities in Carina, as discussed by Humphreys & Kerr (1974). The large range 
of ARq determined within the same data-set reveals just how sensitive the product of these 
constants is to the location in which it is measured. As warned by Lin et al. (1978) one 
should exercise extreme caution when determining ARq from H I or CO terminal velocity 
measurements. 

Recognising the intrinsic errors in determining ARq it is still interesting to compare the 
range of ARq imphed by our data to A measured from proper motions of stars. Comparing 
our values of ARq = 89 — 110 km s"^ with the Hipparchos results for A (Mignard 2000) shows 
that for i?o = 8.5 the lower hmit of our range is consistent with the A values derived for hot 
stars in the disk (e.g. A = 10.9 ± 0.8 km s"^ kpc~^ for stars of type A0-A5 with |6| < 30°), 
but even our upper limit is smaller than the Hipparchos result for giant stars at all latitudes, 
which is A ~ 14.7 ± 1 km s^^ kpc^^. Our upper limit value of A assuming Rq = 8.5 kpc 
is also smaller than the conventionally accepted value of A = 14.5 km s^^ kpc^^ (Binncy 
& Trcmainc 1987). This low value of A derived from the H I terminal velocities compared 
with the locally measured value based on stellar radial velocities and distances supports the 
mounting evidence that the Galactic Center distance is smaller than the 8.5 kpc assumed 
here (Reid 1993; Eisenhauer et al. 2005). 



4.3. Inner Galetxy Rotation Curve 

Using terminal velocities measured from the SGPS data and the circular velocity relation 
that 0(i?) = |fL5K| + ©ol sin/| we have calculated the rotation curve for 3 < i? < 8 kpc. 
Figure 8 shows the rotation curve for the fourth Galactic quadrant derived from SGPS 
terminal velocity measurements. Also plotted on this figure are the rotation curve data 
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points from the first quadrant pubfished by Fich et al. (1989), which were derived from H I 

data presented by Burton & Gordon (1978). We include on Figure 8 two commonly used 
rotation curve fits, the one derived for the outer Galaxy by Brand & Blitz (1993) and the 
one derived for the first quadrant by Burton & Gordon (1978). Note that the Brand & Blitz 
(1993) fit is almost identical to the Fich et al. (1989) rotation curve fit, determined for the 
inner and the outer Galaxy. 

The general trend of the SGPS rotation curve over the full range of plotted radii is 
for increasing rotation velocity with increasing galactocentric radius. In addition, there is 
clearly a great deal of dynamic structure along the locus of subccntral points. The bump 
at i? = 3.4 kpc is attributed to the well-known three kiloparsec arm (van Woerden et al. 
1957). This feature is very well traced by the SGPS terminal velocity data but it is clearly 
not in circular rotation and should not fit by rotation curve models. Other noticeable peaks 
and troughs in the rotation curve have long been attributed to streaming motions associated 
with the spiral arms of the Galaxy (Shane & Bieger-Smitli 1966; Yuan 1969; Burton 1971). 
At these locations the assumption of circular rotation breaks down and the terminal velocity 
curve traces not the overall rotation of the Galaxy, but the dynamical influence of the spiral 
arms. 

From R ^ 3.9 kpc the rotation curve rises rapidly from = 210 km s~^ to 225 km 
at i? ~ 4.4 kpc. This rise is mirrored in the first quadrant rotation curve, also shown in Fig. 
8, although for the first quadrant curve the rise is shifted to larger radii by approximately 
100 pc. At R ^ 5.7 kpc there is a very sharp local minimum of = 211 km s^^ followed by a 
smooth rise towards a maximum of 240 km s~^ near pa 6.7 kpc. There is a clear difference 
between the first and fourth quadrant rotation curves R > 7 kpc. This difference has been 
given as evidence of anomalous velocities associated with streaming motion Sagittarius- 
Carina arm at I ~ 295° (Henderson et al. 1982; Grabelsky et al. 1987). On the other hand, 
asymmetries between the first and fourth quadrant curves have also been explained by Kerr 
(1962) and Blitz & Spergel (1991) as evidence for an outward motion of the Local Standard 
of Rest orbit. If there is an outward motion of the LSR this would apply an offset to the 
derived rotation curve of — 2nocos/, where 11^ is the magnitude of the outward velocity, 
estimated to be 14 km s~^ by Blitz & Spergel (1991). Our current data, because they are 
restricted to the fourth quadrant, are not sufficient to resolve whether there is an outward 
velocity of the LSR so we omit this correction for now. 

The SGPS data, because of their high angular resolution, allow us to clearly resolve 
individual features in the rotation curve. Superposed on the large scale features are a mul- 
titude of small features with velocity amphtudes ~ 2 — 5 km s~^, over length scales of 100 
to 200 pc. One such example lies at ~ 8 kpc and can be attributed to the Carina supershell 
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(Dawson et al. 2006), which has carved a hole out of the disk and disturbed the kinematics 
near the terminal velocity. Presumably some of the other small scale features can be at- 
tributed to discrete objects, but another fraction must represent the random cloud-to-cloud 
motions near the subcentral point. 

4.3.1. A Fit to the Rotation Curve for 3 < R <8 kpc 

As seen in Figure 8, neither of the two commonly used rotation curve formulae are good 
fits to the data. The Brand & Blitz (1993) rotation curve, which is often used for kinematic 
distance estimates, is not a particularly good representation of both the first and fourth 
quadrants. That this curve is a poor fit for the inner Galaxy {R < Rq) is not surprising 
as it was derived to fit the outer Galaxy. When plotted over the inner Galaxy it is clear 
that the Brand & Blitz (1993) fit is much fiatter than the data. Similarly, the Burton & 
Gordon (1978) fit, although a reasonable fit for the first quadrant, is much too steep for 
the fourth quadrant data. Other rotation curves have been derived specifically for the inner 
Galaxy, such as Clemens (1985), but this high order polynomial fit was derived to fit data 
from the locus of subcentral points in the first quadrant. The curve includes fits to the large 
scale velocity deviations known to be associated with spiral arms and is therefore of limited 
applicability to other parts of the inner Galaxy. 

Here we use both the first and fourth quadrant data together to derive a new fit to the 

rotation curve appropriate for 3 kpc < R < 8 kpc at all longitudes. Following Fich et al. 
(1989) and Brand & Blitz (1993) we fit the data with a variety of functional forms, including 
first and second-order polynomials and the power law O/Oq = Ci{R/RoY'^ + <^3> where Ci, 
C2, and C3 are the free parameters to be determined by the fit. The fitting was done using a 
non-linear least-squares algorithm on the combination of the first quadrant and SGPS data, 
including the standard deviations of each. The standard deviation for the first quadrant 
velocities is = 4.5 km s"^ (Fich et al. 1989; Burton & Gordon 1978). The SGPS standard 
deviations were derived from the repeatability of the vlsr determination as discussed in §3.3. 
We let — 1 km s~^ for I < 325° and cr^ = 3 km s~^ for I > 332?5 and ay — 10 km s~^ for 
327?5 < I < 330°. All data were additionally weighted according to sin / to compensate for 
the uneven sampling in this domain. Finally, the first quadrant weights were adjusted by 
the ratio of the number of points in the first and fourth quadrant datasets to compensate for 
the fact that the SGPS data arc much more densely sampled than the first quadrant data. 
We found that the best fit for the range of radii used here, which minimised the absolute 
chi-squared, was a simple linear function, where 



e(i?)/eo = (0.186 ± 0.008) (i?/i?o) + (0.887 ±0.005) 



(6) 
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This fit is shown in Figure 8. We emphasize that this fit is specifically intended for the 
region 3 kpc < R < 8 kpc and is not a good fit to outer Galaxy rotation. This fit, however, 
is generalised for applicability to the inner Galaxy in both the first and fourth Galactic 
quadrants. Although we have assumed Gq — 220 km s~^ in Figure 8, the fit is independent 
of the assumed values for 0o and Rq. 

4-3.2. Velocity Residuals 

From Figure 8 it is clear that there are significant motions that arc not fit by this simple 
rotation curve and the magnitude of these motions should be considered when using this fit 
for kinematic distances. These residuals give an estimate for the magnitude of streaming 
motions in the inner Galaxy. Figure 9 shows the residuals for the SGPS and 1st quadrant 
rotation curves minus the linear fit described above. The peak-to-peak variations in the 
SGPS residuals using the new fit are 10 km s~^ and they have a standard deviation of 5.1 
km s~^. The first quadrant residuals have a larger peak-to-peak variation but this should 
be confirmed by better sampled data. For comparison, the standard deviation of the SGPS 
residuals from the Brand & Blitz (1993) fit are 6.6 km s~^. The magnitude of these residuals 
docs not differ significantly from those estimated by Shane & Bieger-Smith (1966); Burton 
(1971) and others, but is perhaps the most well determined estimation available. 

That the positions of the bumps in the rotation curve relate to spiral structure has 
been long known and STicccssfully demonstrated by applications of density wave theory to 
H I observations of the Milky Way (e.g. Yuan 1969; Burton 1971; Roberts 1969), and other 
galaxies (e.g. Rots 1975; Visser 1980). Indeed, these bumps are the strongest indicators we 
have of the locations of the spiral arms. The problem, however, has been that the linear 
and non-linear versions of the density wave theory predict different relationships between 
the gravitational centre of a spiral arm and the associated velocity peak. Past data from 
the Milky Way has not allowed us to distinguish between these two theories, but with the 
abundance of high resolution H I (McClure- Griffiths et al. 2005; Stil et al. 2006), CO (Luna 
et al. 2006), and infrared data (Benjamin et al. 2003) for the inner Milky Way now available 
it may be time to revisit this topic. This will be the topic of a future paper in this series. 

5. Random Motions Determined Prom The Terminal Velocity Shape 

In addition to providing the offset velocity, Vg, the shape of the terminal velocity profile 
provides information about random gas motions in the ISM. The subcentral point is a useful 
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location for measuring these quantities because velocity crowding forces us to observe a large 
ensemble of H I "clouds" . When we measure the velocity dispersion of H I near the subcentral 
point we measure a mixture of the microscopic (thermal) and macroscopic (the so-called 
cloud-to-cloud velocity dispersion) motions. The terminal profile shape is therefore one of 
the best probes of cloud-to-cloud velocity dispersion in the Galaxy. Our error function fits 
return values for the velocity widths of the error functions, which are related to the velocity 
dispersions of the gas as Avj = a/2(J^. Using the SGPS data beyond the subcentral point 
we have studied the velocity widths both for the dataset as a whole and as a function of 
Galactic radius. 

In Tabic 2 we give values for the fitted parameters for two and three component error 
function fits to the longitude averaged terminal velocity profile. The fitted parameters, a^, 
Vo and Avi, are those defined in Equation 4. The two component fit parameters are given on 
the first line of Table 2 and the three component fit parameters are on the second line. The 
quoted errors are the variation in the parameters that double the reduced of the fit. Note 
that these errors do not account for co-variance of the fitting parameters. As mentioned in 
§3.3 and shown in Figure 3 three components produce a significantly better fit to the data, 
reducing the reduced chi-squared of the fit from over 100 for two components to ~ 4 for 
three components. 

The velocity widths of the two component fit are similar to the Kulkarni & Fich (1985) 
observations of H I emission towards / = 30° and / = 180°, which revealed a narrow com- 
ponent with a velocity dispersion of 0"^, = 6 — 8 km s~^ and a second component with a 
velocity dispersion of 15 — 35 km s~^. Kulkarni & Fich (1985) found that the fits to the 
broad component were not unique and that values between 15 km s~^ and 35 km s~^ fit the 
data equally well. The H I spectra that Kulkarni & Fich (1985) used were from Weaver & 
WiUiams (1973), with an rms noise of ax — 0.38 K, whereas the SGPS longitude averaged 
spectrum probes the tail of the H I profile to ut ~ 50 mK. With the better sensitivity of 
the longitude averaged SGPS spectrum we can clearly see that the data are better described 
with three velocity components. Compared to the two component fit, the two narrowest 
components of the three component fit arc smaller. We sec that in the three component 
fit the amplitudes of the first two components are similar and the amplitude of the third 
component is ~ 10% of the first. 

That the terminal velocity profile requires three components is perhaps not surprising. 
H I observed in emission is a mixture of warm (T ~ 6000 K) and cool (T = 50 — 100 K) 
kinetic temperatures. If we ignore the effects of the velocity gradient, the observed velocity 
dispersion, a„, represents a combination of the thermal velocity width, ath, and the random 
velocity of the cloud ensemble (cloud-to-cloud motions), (7cc, such that the total velocity 
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width is Av = V^a^ = \/H^^h+~^- ^ 

we assume thermal velocity dispersions for the 
cold and warm gas of ~ 1 km and ~ 7 km s~^, respectively, then the widths of the 
two narrow components give cloud-to-cloud velocity dispersions of ~ 4.4 km and ~ 5.2 
km s~^, which are identical to within the errors. This suggests that the cold and warm H I 
share the same bulk motions. 



5.1. Velocity Dispersion Variations with Galactocentric Radius 

Kulkarni & Fich (1985) suggested that the wide velocity dispersion might be related 
to shocked clouds. For this reason they chose to examine a profile at / = 30° where the 
line-of-sight runs tangent to the R = 5 kpc annulus where supernova remnants show a peak. 
Additionally, Koo and collaborators (i.e. Koo & Kang 2004) have suggested that extended 
forbidden-velocity line wings may be associated with individual supernova remnants, whose 
shocks drive H I to high velocities. Given these suggestions it is interesting to explore whether 
there may be a significant longitude dependence to the velocity dispersion. 

Using the full range of the SGPS data wc explore whether the velocity widths of the 
three error functions change with position in the Galaxy. For this analysis we include data 
in the range 285° < I < 340°, which excludes data from the inner 3 kpc of the Galaxy 
where the ISM environment is different from the rest of the disk. We have also excluded 
data for I < 285° because at these longitudes the velocity gradient, dv/dr, is near zero for a 
particularly long path length, leading to saturation in the H I profiles. Over the longitude 
range 285° < / < 340° there are ~ 1520 independent samples. However, to explore the 
broad velocity component we are interested in brightness temperatures extending as low as 
~ 250 mK, so we must average the data into longitude bins of ~ 1?4 to achieve sufficient 
signal-to-noise. Note that this analysis necessarily precludes us from finding possible missing 
supernova remnants by their forbidden-velocity line-wings because we average over large 
spatial ranges. For each longitude bin we have fit the terminal velocity curve with a sum of 
three error functions, as defined in Equation 4. The difference between three and two error 
functions is important here as a fit with only two components produces a consistently wider 
first component than does the three component fit. 

In Figure 10 we plot the widths of the three error functions, A.vi{R), Av2{R), and 
Av^{R) as a functions of galactocentric radius. The error bars are estimated as the variation 
in the parameter that results in a doubling of the reduced chi-squared. In this plot the three 
different symbols correspond to the three different components. The mean values of the 
three widths are: 6.5 ±2.1 km s^^, 12.8 ±3.2 km s^^ and 25.5 ±7.9 km s^^, which agree well 
with the values found for the fit to the data averaged over the entire longitude range. We 
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have fit both of the narrowest components as functions of radius, shown in Figure 10. The 
slope of the narrowest component is consistent with zero: Avi{R) = (6 ± 1) + (0.2 ± 0.2)R 
km s~^. There is only a very slight slope to the second component, which goes as Av2{R) = 
(13±1) — (0.3±0.2)i? km s~^. The widest component does not show a linear trend but shows 
a significant increase at small Galactic radii. This increase is accompanied with increasing 
errors because at large longitudes the tail of the H I emission extends to the edge of the 
SGPS bandwidth, making it difficult to fit the long tail of the H I profile. Despite the larger 
errors, the increase appears real. 

Burton (1971) also examined the velocity dispersion versus Galactic radius, fitting his 
single component velocity dispersion with a function of the form a{K) = 9 — 0.47? km s^^. 
Our fit for the narrowest component is consistent with no gradient in the velocity dispersion 
across the disk and the middle component shows a shallower slope than Burton's result. At 
first our result appears inconsistent with observations of velocity dispersion across the disk 
of face-on external galaxies, which, like Burton (1971) show that inside the stellar disk the 
narrow component velocity dispersion increases with decreasing radius (Dickey et al. 1990). 
The resolution of this apparent conflict, however, may lie in the fact that both studies were 
based on single component fits, which will produce wide velocity dispersions if the profiles 
would be better fit by two or three Gaussian components. From our work it seems that the 
narrow components are effectively constant with radius while the broad component increases 
with decreasing radius. 

In addition to the large-scale trends discussed above, all three components show some 
small-scale features. At i? ~ 5 kpc the wide component increases by a factor of ~ 2 over 
a width of ~ 500 pc. This increase is also refiected in increases in both of the narrow 
components. We see that at it! ~ 6.3 kpc and R ~ 6.8 kpc the fit effectively reduces to a 
two component model. It may be noteworthy that the narrowest two components track each 
other very well with Galactic radius, with all fiuctuations observed in both components. This 
observation supports our suggestion in the previous section that the warm and cold H I share 
the same bulk motions, which may be either large-scale streaming motions or cloud-to-cloud 
motions. The separation between the two widths is simply the difference in the thermal 
motions, not the macroscopic motions which appear to effect warm and cold gas equally. 

Some observations of face-on galaxies suggest that there may be some correlation be- 
tween velocity dispersion and spiral structure (e.g. Braun 1997). Similarly, the original 
supposition of Kulkarni & Fich (1985) that led them to look for "fast" clouds at / = 30° 
was that there may be correlation between fast clouds and areas rich with supernovae, such 
as spiral arms. At first glance Figure 10 seems to support this latter supposition because 
all velocity components increase near R — 5 kpc. It is tempting to look deeper into our 
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fitted velocity dispersions for correlation with spiral structure but Burton (1971) reminds 
us that at the subcentral point the measured velocity dispersion can be related to the path 
length contributing to the velocities near the subcentral point, in other words the inverse 
of the velocity gradient, \dv/dr\. In the presence of streaming motions \dv/dr\ is altered 
from the value predicted by circular rotation. An increase in velocity dispersion is therefore 
not necessarily indicative of an increase in the cloud-to-cloud dispersion, but could simply 
indicate an increase in the path length contributing to the emission near the terminal veloc- 
ity. Burton showed that for linear density wave theory streaming motions should increase 
the path length near the terminal velocity just outside the spiral arms and suggested that 
there was an indicative increase in his measured dispersions. However, the Burton (1971) 
analysis was hampered by single Gaussian fits, which, as he pointed out, were extremely 
difficult at longitudes where the velocity component became so broad with a flat top that it 
was better flt by two separate components. Fitting error functions to the tail of the emission 
beyond the terminal velocity provides some protection against this problem but there is no 
doubt that in the presence of streaming motions the measured velocity dispersion will change 
from that measured given purely circular rotation. Unfortunately, once again the location of 
increases relative to spiral structure depends on the density wave theory applied, linear or 
non-linear. Without knowing the gravitational positions of the spiral arms it is impossible 
to say whether the increases in velocity width that we observe in all components at i? ~ 5 
kpc and in the largest component at < 4 kpc are indicative of a change in dv/dr or a 
change in the cloud-to-cloud velocity dispersion. 

The velocity gradient also plays a role in the mixture of gas-phase components measured 
beyond the terminal velocity. In the simplest case this can be understood if we consider that 
for emission beyond the subcentral point and in a given velocity interval, 6v, the path length 
contributing to the emission is 6r = 5v{dv/dr)^^. We can therefore see that for gas with 
a narrow velocity dispersion only gas very near to the subcentral point contributes to the 
emission beyond the terminal velocity, whereas gas with a large velocity dispersion can 
be distributed over a large path length near the subcentral point and still contribute to the 
emission beyond the terminal velocity. Precisely at the subcentral point the velocity gradient, 
\dv/dr\ — so that it is the second derivative of velocity with distance, d'^v/dr'^, that 
determines the distance along the line of sight that contributes to an interval of velocity near 

the subcentral point os 5r — ^J2 5v \d'^v / dr'^\~^ , i.e. the second term of a Taylor expansion 
of v{r) for r ~ tq. The larger the second derivative, the smaller the line of sight contributing 
to a given spectral interval. Even in pure circular rotation the second derivative of velocity 
near the subcentral point changes as a function of longitude so we expect that the proportion 
of phases represented beyond the terminal velocity will change with longitude. As discussed 
above, the velocity gradient is also sensitive to streaming motions so understanding the 
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mixture of phases with radius is an extremely comphcated issue. We defer that discussion 
to a future paper, which will explore the spiral structure more fully. 

6. Conclusions 

For the A21-cm pioneers who first measured the Milky Way rotation curve 50 years 
ago, the objective was to determine the radial variation of the gravitational potential, and 
so to measure the distribution of gravitational mass. In this they were quite successful 
(Kwcc ct al. 1954; Schmidt 1956) and only relatively minor adjustments have been made 
since by breaking the Galaxy down into its disk and spheroidal components (Caldwell & 
Ostriker 1981). Refinements such as those discussed in this paper do not alter the overall 
result for the mass of the Galaxy. The significance of the high angular resolution and high 
precision measurement techniques developed here have more to do with departures from 
circular rotation than with the overall mass versus radius function. We have fit a new 
rotation curve to the region 3 < < 8 kpc using our data from the fourth quadrant 
together with a previously published rotation curve for the first quadrant. This generalised 
fit attempts to overcome problems related to quadrant specific fits and is applicable for 
all of the inner Galaxy. Using this fit we determine the velocity residuals for the fourth 
quadrant confirming with greater precision that the streaming motions have amplitudes of 
~ 10 km s^^. Although the wiggles evident on Figure 9 have long been interpreted in terms 
of spiral structure their precise location with respect to the gravitational centre of spiral 
arms is theory dependent. We have future work planned to combine high precision H I and 
CO surveys with observations of star-forming regions to better quantify the spiral structure 
in the inner Galaxy. 

In addition to large scale variations in the magnitude and direction of the streamlines 
of the interstellar gas as it rotates around the Galactic Center, some of the small scale 
fiuctuations measured here may be due to non-gravitational perturbations that change the 
shape of the line profile beyond the terminal velocity. In §3 and §5 we carefully discussed the 
velocity offset derived from error function fitting to the profile shape beyond the terminal 
velocity, because the random motions reflected in this shape must be separated from large 
scale variations in the rotation velocity. Thus in order to trust the interpretation of the 
wiggles on Figure 9 as being due to spiral arms, we must first eliminate wiggles in Vg similar 
to those in Avi, Av2, and Av^ on Figure 10. 

We have presented a successful technique for measuring the terminal velocity of H I 
in the inner Galaxy. As shown in Figures 4 and 5, the technique of fitting complementary 
error functions to the H I terminal velocity is quite robust, giving results that are consistent 
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to less than the width of one velocity channel over much of the longitude range of the 
survey. Our technique is effective at overcoming the biases that might be introduced by 
determining the terminal velocity using brightness temperature threshold methods. Our fits 
show that H I terminal velocities estimated based on a brightness temperature threshold 
without any correction will over-estimate the terminal velocity by 5 — 10 km s~^. Although 
the CO velocity dispersion is smaller than for H I, terminal velocities estimated from CO 
emission using a threshold technique are also not immune from over-estimation by amounts 
comparable to molecular cloud velocity dispersion. Comparison of the H I and Luna et al. 
(2006) ^^CO terminal velocities has shown that the CO velocities have more scatter than 
those presented here from H I. We attribute some of this difference to the different volume 
filling factors of H I and CO and suggest that H I measurements of the terminal velocity 
may be more robust than CO because they detect a larger ensemble of clouds. 

The technique of fitting the H I terminal velocity profile has the added benefit that it 
returns information about the velocity width of the H I beyond the terminal velocity. The 
fitted velocity width combines information about the thermal width of the gas, the inherent 
velocity dispersion as caused by cloud-to-cloud motions, and changes in the velocity gradient, 
which may be caused by streaming motions associated with spiral arms. Unravelling these 
effects at specific positions within the inner Galaxy is complicated and something we hope to 
address in a future paper. Averaging over all longitudes to mitigate against local variations 
in the velocity gradient we have shown that the H I terminal velocity shape is best fit by the 
combination of three error functions, with widths of 6 km s~^, 12 km s~^, and 26 km s~^. 
Examining the three velocity widths as a function of Galactic radius we found that the 
widths of the narrower two components track each other very well at all Galactic radii and 
show no appreciable trend with radius. We suggest that the agreement between the narrow 
components is evidence that the two thermal phases share the same bulk motions as well 
as variations in the velocity gradient. Unlike the two narrow components, the third, widest, 
component appears to increase with decreasing Galactic radius for R < 4 kpc. It is unclear 
at this time whether this increase refiects a true increase in the cloud-to-cloud motions or 
streaming motions. 
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Fig. 1. — An illustration of the effect of continuum absorption on the H I spectra. The left 
hand panel shows a grey scale map of continuum sources near longitude / = 305?25. The 
vertical line spanning +0?5 > b > — 0?5 shows the area covered in a single latitude average. 
The latitude-velocity plot in the center panel shows the H I line intensity along this line, 
with velocity as the horizontal axis and grey scale indicating brightness temperature. The 
white areas have Tb < because the continuum has been subtracted from the spectral line 
cubes. The right hand panels show two spectra, at 6 = +0?2 and b — — 0?2, shced from the 
latitude velocity plot. The upper one is polluted by absorption toward the continuum and 
must be excluded from the calculations. 
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Fig. 2. — Distribution function of the brightness temperature measured ai v — —61.83 
km / = 305?25, for latitudes +0?5 >b> -0?5, which is the range shown on Figure 1. 
The velocity displayed here is the threshold crossing velocity for this longitude, meaning it 
is the most negative velocity for which the mean of the brightness temperature distribution 
is above the 10 K threshold. The mean, median, 25th percentile, and 75th percentile of the 
distribution are indicated by the marks labelled 1, 2, 3, and 4 respectively. The shaded pixels 
are eliminated from the distribution before computing the statistics because they correspond 
to positions where the continuum is above the cutoff of 0.5 K. Distribution functions like this 
are analysed at every longitude, and the four statistics are used to compute four alternative 
velocity shifts, Vt{l), where each one crosses the 10 K threshold. 
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Fig. 3. — Error function fits to the longitude averaged terminal velocity profile. Panel a 
shows a two component fit to the profile while panel h shows a three component fit to the 
profile. Clearly the three component fit is much better able to account for the tail of the 
profile. 
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Fig. 4. — The effect of correcting tlie tlireshold crossing velocity by tlie velocity offset 
from the error function fits. The four crosses at each longitude are the average values of Vt in 
each longitude interval (2?5), for the four statistics (mean, median, 25th and 75th percentile) 
used to find the threshold crossings. The 25th percentile gives consistently lower negative 
velocity values, while the 75th percentile gives higher values, while the mean and median 
give values of Vt that are generally quite close together. But after adding the corresponding 
Vo all the values agree much better, as indicated by the curves that join values of vlsr for 
each longitude bin. The dotted curve is for the 75th percentile, the dashed curve for the 
25th percentile, and the two solid curves are for the mean and median. All four agree to 
better than 0.5 km s~^ for longitudes less than about / = 325°, as discussed in the text. 
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Fig. 5. — The difTcrcncc between choosing thresholds of 10 K and 20 K in brightness tem- 
perature. As on Figure 4 the crosses show the threshold crossing velocities, this time for the 
two thresholds both using the mean value (average over latitude) as the statistic. The higher 
negative velocities correspond to the lower threshold, as shown on fig. 4. After applying the 
correction Vg found from erfc fitting, the two thresholds give the same result to less than 0.5 
km s~^ for longitudes less than 325°, as discussed in the text. This demonstrates that the 
final function VLsniP) is independent of the arbitrary choice of threshold over most of the 
longitude range of the survey. 
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Fig. 6. — H I terminal velocities, vlsr, and the offset velocity, Vo, versus Galactic longitude 
for the fourth quadrant of the Milky Way. 
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Fig. 7. — H I and CO terminal velocities versus | sin/| for the fourth quadrant of the Milky 
Way. CO terminal velocities are from Luna et al. (2006) [L06]. The H I terminal velocities 
have been fit with a hnear function of the form vlsr = ©od sin/| — 1) + v'. The slope of 
the fit is dependent on the range of sinl used in the fitting. The dashed fine shows a fit 
over the range 0.35 < | sinZ| < 0.95 and the sohd fine shows the fit over the restricted range 
0.8 < I sinZ| < 0.95. The fitted values are given in the legend. 
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Fig. 8. — Rotation curve for the fourth quadrant of the Milky Way interior to the Solar circle. 
The crosses are the SGPS data points and the open circles are the first quadrant values from 
Fich et al. (1989). Error bars are not shown on this plot to reduce confusion. Estimates of the 
errors are included in the fitting and are 4.6 for I > 0° (Fich et al. 1989), 1 km s~^ for I < 325°, 
3 km s-^ for I > 332?5 and 10 km s"^ for 327?5 < I < 332°, as described in the text. The 
dotted line is the commonly used Brand & Blitz (1993) rotation curve derived for the outer 
Galaxy. The dashed-dotted line is the Burton & Gordon (1978) polynomial fit to the first 
quadrant rotation curve. This fit has been scaled from the original assumption of Rq = 10 kpc 
and Oq = 250 km s~^ to Rq = 8.5 kpc and ©q = 230 km s~^. We found it necessary to use 
Go = 230 km s~^ rather than 220 km s~^ to place the peak near the observe peak velocities. 
The sohd fine is the best fit to the ensemble of the first and fourth quadrant (SGPS) data 
over the range 3 < < 8 kpc given by e(i?)/eo = 0.186 ± 0.008(i?/i?o) + 0.887 ± 0.005, 
where ©o — 220 km and Rq — 8.5 kpc. 
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Fig. 9. — Residuals from the rotation curve fit presented in Figure 8. The standard deviation 
for the SGPS residuals is 4.11 km s~^, whereas the standard deviation for residuals from the 
Brand & Blitz (1993) curve was 6.6 km s~^. 
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Fig. 10. — Widths of the three velocity components fitted to the terminal velocity profile 
as a function of Galactic longitude. These data are binned by 1?4 of Galactic longitude. 
The plotting symbols for the narrow, middle and wide components are circles, triangles and 
boxes, respectively. The error bars are the amount of variation in a parameter required to 
double the chi-squared of the fit. The narrowest components are fit with linear functions 
of the form Avi{R) = (6 ± 1) + (0.2 ± 0.2)R km s'^and Av2{R) = (13 ± 1) - (0.3 ± Q.2)R 
km s~^. 
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Table 1. Measured H I terminal velocities, vlsr, as a function of Galactic longitude, /, 
from the SGPS. The third column gives the terminal velocities, Vmod, predicted by the 
rotation curve defined in Equation 6. Note- the full version of this table appears in the 
electronic edition of the Journal, an abbreviated form is presented here. 



I (deg) 


VLSR (km s 1) 


Vmod (km s 1) 


339.695 


-135.22 


-133.00 


339.628 


-135.14 


-132.80 


339.562 


-135.44 


-132.61 


339.495 


-136.17 


-132.41 


339.428 


-135.59 


-132.22 


339.362 


-137.87 


-132.02 


339.295 


-136.53 


-131.83 


339.228 


-135.63 


-131.63 


339.162 


-135.47 


-131.44 


339.095 


-134.11 


-131.24 



Table 2. Fitted values for velocity widths of two and three component error function fits 
to the terminal velocity curve averaged over the longitude range 272° < I < 340°. 



i;o(kms-i) ai (K) Ai;i (km s"!) a2 (K) Ad2 (km s^i) as (K) At;3 (km s'^) 



2- component fit 7.7 ±0.3 45.4 ±3.6 8.6 ± 0.3 8.1 ±0.8 21.2 ± 1.3 

3- component fit 7.0 ±0.1 28.6 ±1.1 6.3 ±0.1 19.5 ±0.4 12.3 ± 0.3 3.1 ±0.1 25.9 ± 0.5 



